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Abstract 

Neutrino oscillations are well established and the relevant parameters determined with good 
precision, except for the CP phase, in terms of a unitary lepton mixing matrix. Seesaw extensions of 
the Standard Model predict unitarity deviations due to the admixture of heavy isosinglet neutrinos. 
We provide a complete description of the unitarity and universality deviations in the light neutrino 
sector. Neutrino oscillation experiments involving electron or muon neutrinos and anti-neutrinos 
are fully described in terms of just three new real parameters and a new CP phase, in addition to the 
ones describing oscillations with unitary mixing. Using this formalism we describe the implications 
of non-unitarity for neutrino oscillations and summarize the model-independent constraints on 
heavy neutrino couplings that arise from current experiments. 
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I. INTRODUCTION 


Neutrino masses, without which current neutrino oscillation data can not be under¬ 
stood [1], are here to stay [2]. It has been long noted that small neutrino masses can 
arise from an effective lepton number violation dimension-five operator O 5 oc , which 

may arise from unknown physics beyond that of the SU(3)c ® SU(2)l ® U(1)y model. Here 
L denotes one of the three lepton doublets and d* is the standard model scalar doublet [3]. 
After electroweak symmetry breaking takes place through the nonzero vacuum expectation 
value (vev) (d)) such operator leads to Majorana neutrino masses. In contrast to the charged 
fermion masses, which arise directly from the coupling of the scalar Higgs, neutrino masses 
appear in second order in (d)) and imply lepton number violation by two units (AL = 2) at 
some large scale. This fact accounts for the smallness of neutrino masses relative to those 
of the standard model charged fermions. This is all we can say from first principles about 
the operator O 5 in Fig. [Tj In general we have no clue on the mechanism giving rise to O 5 , 
nor its associated mass scale, nor the possible details of its flavour structure. 


L 



FIG. 1: Dimension five operator responsible for neutrino mass. 



S[/(3) ® 5C/(2) ® [/(l) 

La — (^a? la) 

(1,2,-1/2) 

(1,1,1) 

Qa — (Ua, da) 

< 

da 

(3, 2,1/6) 

(3,1,-2/3) 

(3,1,1/3) 


(1,2,1/2) 


TABLE I: Matter and scalar multiplets of the Standard Model. 

One may assume that O 5 is induced at the tree level by the exchange of heavy “messenger” 
particles, whose mass lies at a scale associated to the violation of the global lepton number 
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symmetry by new physics, beyond that of the SU(3)c ® SU(2)l ® U(1)y model, 


m„ = Xq 


{31 

Mx' 


where Aq is some unknown dimensionless constant. For example gravity, which in a sense 
’’belongs” to the SM, could induce the dimension-hve lepton number violation operator 
(^5 ais]. In such a minimalistic scenario [6] the large scale Mx in the denominator is the 
Planck scale and hence the neutrino mass that results is too small to account for current 
neutrino oscillation data. Hence we need genuine “new physics” in order to generate neutrino 
masses this way. 

Neutral heavy leptons (NHL) arise naturally in several extensions of the Standard Model. 
Their possible role as messengers of neutrino mass generation constitutes one of their 
strongest motivations and a key ingredient of the type-I seesaw mechanism [ZHn] in any of 
its variants. If realized at the Fermi scale [T21 - I20] . it is likely that the “seesaw messengers” 
responsible for inducing neutrino masses would lead to a variety of phenomenological impli¬ 
cations. These depend on the assumed gauge structure. Here for definiteness and simplicity, 
we take the minimal SU(3)c ® SU(2)l ® U(1)y structure which is well tested experimentally. 
In this case one can have, for example. 


1. Light isosinglet leptons are usually called “sterile”. If they lie in the eV range they 
could help accommodate current neutrino oscillation anomalies |2IIE2] by taking part 
in the oscillations. Sterile neutrinos at or above the keV range might show as distor¬ 
tions in weak decay spectra [23] and be relevant for cosmology 121 . 

2. Heavy isosinglet leptons below the Z mass could have been seen at LEP I p31f?F] . 
Likewise, TeV NHLs might be seen in the current LHC experiment, though in the 
latter case rates are not expected to be large in the SU(3)c ® SU(2)l 0 U(1)y theory. 

3. Whenever NHL are too heavy to be emitted in weak decay, processes, the correspond¬ 
ing decay rates would decrease, leading to universality violation [2H|- 

4. The admixture of NHL in the charged current weak interaction would affect neutrino 
oscillations, since they would not take part in oscillations. These would be effectively 
described by a non-unitary mixing matrix 123. 

5. If Majorana-type, NHL would modify rates for lepton number violation processes such 
as neutrinoless double beta {Oh'(3f3) decays through long-range (mass mechanism), as 
well as induce short-range contributions [30f[32] . 

6. NHL would induce charged lepton flavour violation processes [29l |33|. However the 
corresponding restrictions depend on very model-dependent rates. 
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In what follows we consider the generic structure of the lepton mixing matrix relaxing 
the unitarity approximation We show that their most general form is factorizable, so that 
current experiments involving only electron and muon neutrinos or anti-neutrinos can be 
effectively described in terms of just three new real parameters and one new CP violation 
phase. We illustrate how these parameters affect oscillations and discuss the main restric¬ 
tions on such generalized mixing structure that follow from universality tests. For logical 
completeness we also present a brief compilation of various model-independent constraints 
on NHL mixing parameters within the same parametrization, including those that follow 
from the possibility of direct NHL production at high energy accelerator experiments. 


II. THE FORMALISM 


Isosinglet neutral heavy leptons couple in the weak charged current through mixing with 
the standard isodoublet neutrinos. The most general structure of this mixing matrix has 
been given in the symmetric parametrization in Ref. [8]. Here we consider an equivalent 
presentation of the lepton mixing matrix which manifestly factorizes the parameters asso¬ 
ciated to the heavy leptons from those describing oscillations of the light neutrinos within 
the unitarity approximation. Here we present its main features, details are given in the 
appendix 

For the case of three light neutrinos and n — 3 neutral heavy leptons, one can break up 
the matrix describing the diagonalization of the neutral mass matrix as [35] 


U 


nxn 



( 1 ) 


where iV is a 3 x 3 matrix in the light neutrino sector, while S describes the coupling param¬ 
eters of the extra isosinglet states, expected to be heavy (for a perturbative expansion for 

jjnxn ggg 

M)- As shown in the appendix, the matrix N can be expressed most conveniently ^ 


as 



f Oill 

0 

0 \ 

N = N^^U = 

Ct21 

“22 

0 


V “31 

“32 

“33 / 


( 2 ) 


^ In sections II-VI we mainly consider isosinglet neutrinos above 100 GeV or so, hence too heavy to take 
part in oscillations or low energy weak decay processes. 

^ We consider stable neutrinos, neutrino decays were discussed, for instance, in Ref. [34| . 

^ There are other forms for the light-neutrino mixing matrix, where the pre-factor off-diagonal zeroes 
are located at different entries. However Eq. ([^ is the most convenient to describe current neutrino 
experiments. 
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where U is the usual unitary form of the 3x3 leptonic mixing matrix probed in neutrino 
oscillation studies ^ corrected by the left triangle pre-factor matrix, characterizing 

unitarity violation. 

Note that Eq. (|^ provides a most convenient, general and complete description of the 
propagation of solar, atmospheric and terrestrial neutrinos from reactors, radioactive sources 
and accelerators beams, relaxing the unitarity approximation. Due to the zeroes in the first 
two rows of the pre-factor matrix in Eq. ([^ it is clear that the only extra parameters beyond 
those characterizing unitary mixing are four: the two real parameters an and 022 plus the 
complex parameter 021 which contains a single CP phase. Indeed the existence [3Z] and 
possible effects [38] extra CP phases associated to the admixture of NHL in the charged 
leptonic weak interaction had already been noted in the early paper in [8|. The new point 
here is that, despite the proliferation of phase parameters, only one combination enters 
the “relevant” neutrino oscillation experiments. This holds irrespective of the number of 
extra heavy isosinglet neutrino states present. Other studies, such as [39IHT] . appear as 
particular cases with a hxed number of extra heavy isosinglet neutrino states, any of which 
can be expressed in terms of the same set of parameters aij. Similarly, the matrix U may 
be expressed in different ways, such as in PDG form or in our fully symmetric description, 
particularly useful for phenomenological analyses. The diagonal elements, an, are real and 
expressed in a simple way as 


ail 

^In ^ln—l^ln—2 • 

• • Cl4, 

OC 22 

C2n C2n—1^2n—2 • 

• • C 24 , 

tt33 

C372 C^rL—1^3n—2 • 

• • C 34 , 


in terms of the cosines of the mixing parameters [S], Cij = cos 6ij. 

Now the off-diagonal terms a 2 i and a 32 are expressed as a sum of n — 3 terms 


(3) 


0^21 — C 2 n C 2 n -1 ■ ■ ■ C25 ^ 24^14 + C 2 n • • • C26 ^ 25^15 Ci 4 + • • • + ^ 2 n^ln Dn-l Cin -2 ■ ■ ■ C14 , 

0^32 = C 3 n C 3 n -1 ■ ■ ■ C35 rj 34 fj 2 A + C3„ . . . C36 ^ 735^25 C24 + • • • + ^ 3 n^ 2 n C 2 n -1 C 2 n -2 ■ ■ ■ C 24 ,( 4 ) 

where rjij = sinOij and its conjugate fjij = —sinOij contain all of the CP violating 

phases. Finally, by neglecting quartic terms in sin^jj, with j = 4 , 5 , • • • one finds a similar 
expression for asi, 

•^31 = C3„ C3„_1 . . . C35 7734C2 4?7i 4 -f- €3^ ■ ■ ■ C3 q r]33C2 ^fii^ €14^ + ... 

T ^ 3 n^ 2 n^ln ^ln—1 C-ln—2 • • • C44 . ( 5 ) 


^ As discussed in Ref. |dfi| . this may, for example, be parameterized in the original symmetric way or 
equivalently as prescribed in the Particle Data Group. 
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In summary, by choosing a convenient ordering for the products of the complex rotation 
matrices Uij (see appendix), one obtains a parametrization that separates all the information 
relative to the additional leptons in a simple and compact form, containing three zeroes. We 
will now concentrate on this specihc parametrization. 


III. NON-UNITARY NEUTRINO MIXING MATRIX 


Given the above considerations and the chiral nature of the SU(3)c ® SU( 2 )l ® U( 1 )y 
model, we notice that the couplings of the n neutrino states in the charged current weak 
interaction can be described by a rectangular matrix |8] 

K=(^N S), (6) 

with iV a 3 X 3 matrix described by Eq. ([^ and S' a 3 x (n — 3) matrix. This can be 
parametrized in the symmetric form or as prescribed in the Particle Data Group. The 
relative pros and cons of the two presentations are considered in Ref. [36] . 

The presence of extra heavy fermions that mix with the active light neutrinos would 
imply the effective non-unitarity of the 3x3 light neutrino mixing matrix, hence modifying 
several SM observables. For example, note that the unitarity condition will take the form 


KK^ = NN^ + SS^ = J, 


(7) 


with 


/ a 


NN^ = 


2 

11 


0^220^32 + CI 21 CI 31 


\ 

'V 


( 8 ) 


^ 11^21 

OillO '21 Q!22+|ci21P 

V CI11CI31 CI22CI32 + CtSlQ^^l CI33 + |ct 3 i|^ + ICI32I 

We will show that, with the parametrization discussed here, one can, at least in principle, 
introduce all of the information of the extra n — 3 states into the aij parameters in a simple 
compact form. The method is completely general and includes all the relevant GP phases. 
In what follows we will consider different direct or indirect tests of the existence of the extra 
heavy fermions, expressing the relevant observables in terms of these parameters, in order 
derive the relevant constraints. 


IV. UNIVERSALITY CONSTRAINTS 

First one notes that if, as generally expected due to their gauge singlet nature, the 
heavy leptons can not be kinematically emitted in various weak processes such as muon or 
beta decays, these decays will be characterized by different effective Fermi constants, hence 
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breaking universality. One can now apply the above formalism in order to describe the 
various weak processes and to derive the corresponding experimental sensitivities. We hrst 
discuss the universality constraint, already reported in the literature |28[ in order to 

cast it within the above formalism. Comparing muon and beta decays one hnds 


and 


Gm = Gf ViNN^)iiiNm)22 = Gf + 

Gp = Gf = Gf 


«ll(«22 + |a2iP), 

(9) 

\[^V 

(10) 


Therefore, all the observables related to Fermi constant will be affected by this change, for 
instance, the quark CKM matrix elements [12]. In particular, the CKM matrix elements Wd 
and Vus are proportional to G^. These matrix elements are measured in /3-decay, iCes decay, 
and hyperon decays. The effect on G^, therefore, modihes Wi and the unitarity constraint 
for the hrst row of the CKM is now expressed as 


El'"" 


2=1 


G^ 

G,, 


GFViNm)u 


GF^/{Nm),,{Nm)^2j iNN%2 


( 11 ) 


where the Eq. (|^ has been used in the last equality. Following the previous equation one 
gets [SO] : 

1 


El'"" 


2=1 


0^22 + \0l2l \ 


= 0.9999 ± 0.0006, 


( 12 ) 


and, therefore, 1 — (iViVl )22 = (*S'S'1)22 = 1 — <^22 “ |Q' 2 iP < 0.0005 at la. 

There are other universality tests that give constraints on these a parameters. For ex¬ 
ample, universality implies that the couplings of the leptons to the gauge bosons are havor 
independent, a feature that emerges in the the standard model without heavy leptons. In 
the presence of heavy isosinglets, these couplings will be havor dependent; the ratios of these 
couplings can be extracted from weak decays and they are expressed as 


{NN% 


a = 1, 3 . 


(13) 


{NNi)22 

For a = 1, this ratio can be constrained by comparing the experimental measurement 
and the theoretical prediction of the pion decay branching ratio m- 

F(7r^ 


Rtt = 




F(7r+ -)■ fi+u)' 


(14) 


One obtains @51 El]: 

R. {NN^)n 


a 


= 


11 


(1.230 ±0.004) X 10 


-4 


RSM {Nm )22 ai2±l«2iP (1.2354 ± 0.0002) X 10-4 


= 0.9956 ± 0.0040 (15) 
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FIG. 2: Constraints on the deviations from unitarity. 


which implies 1 — ah < 0.0084 at la for the least conservative case of ah + |tt 2 ip = 1. 
This procedure was adopted in Ref. |17j. However, in general [(iViVl) 22 ] 7 ^ 1, and it can be 
estimated using the unitarity constraints on the CKM matrix discussed above. Combining 
both constraints (from Eqs. (12) and (15)) we obtain the results shown in Figure]^ restricting 
the parameter combinations shown in the plot. These translate in the constraints 


1-ah < 0.0130, 

1- ah - |a2iP < 0.0012 , 


(16) 


at 90% C.L. for 2 d.o.f. One can make use of a third observable in order to have constraints 
for every independent parameter. This will be discussed in the next section. 

For the sake of completeness we now show the constraints coming from the fi — r univer¬ 
sality which, using Fq. (13), give the bound: 


(A^Alt)33 

{NNh22 


0.9850 ±0.0057. 


(17) 


This implies 1 — < 0.0207 at la for the least conservative case of 

{38^22 = 0. The experimental value was taken from Ref. [52]. We now turn to neutrino 
oscillations. 








V. NON-UNITARITY EFFECT ON NEUTRINO OSCILLATIONS 


In this section we focus on neutrino oscillation experiments. First we obtain general 
expressions for neutrino survival and conversion probabilities in this parametrization and 
confront them with the existing experimental data. The general expressions will be relatively 
simple, especially if we neglect cubic products of a 2 i, sin 6 ^ 13 , and sin(^^), which is a 
reasonable approximation for many applications. The results of this approach for the three 


probabilities discussed in this section are shown in Eqs. (21), (27) and (33). 

For the case of the muon neutrino conversion probability into electron neutrino we have: 


p = 

/ie 






j>i 

3 

j>i 




sm 


AE 


Arrij^L 


2E 


(18) 


And now, instead of the usual unitarity condition for the 3x3 case, we must use the condition 
given in Eqs. ([^ and (|^, arriving to the expression 



Using Eq. (|^ one can substitute the values of Nai in terms of Uai and aij to obtain 

Pfie = aiilaalP 1 1 - 4^ |[/ej 

V j>i 

3 

- (aiia22)'4 5];i?e[UTf/, 

j>i 
3 

+ (0110:22)^2 sin 

j>i 
3 

- 40:^0:22 Re [a2i\Uei\^UljUej + sin^ 

j>i 

3 f Arn?-L\ 

+ 2 al^a 22 ^ Im [a 2 i\Uei\^U*^jUej + a* 2 i\Rej?U^,iUl^ sin ( j . ( 20 ) 

j>i h / 


/A^\ 


/A^\ 



Substituting the terms Uai in our parametrization, and neglecting cubic products of 0 : 21 , 
sin 013 , and Am|i, one obtains 
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P^le — {<yua22Y + «iia22 |021 |-P;[e + 0^1l\o^2l\'^ , 
where we have denoted the standard three-neutrino conversion probability Pj^g^ as 


( 21 ) 


p3x3 _ A 

{le ^ 


, 2/1 a -2 a „• 2 


cos 6*12 COS 6*23 sin 6*12 sin 






2/3 „• 2/3 „• 2 /3 ■„2 f 


-|- cos^ 6*13 sin^ 6*13 sin^ 6*23 sin 


V 


( 22 ) 


T j 


n,±j \ . I Am^^L\ f Arrin-fL 

' sin —cos --41- 1 


4B. 


V 


\ 4E, 


123 , 


while P^g refers to a term that depends on the 3x3 mixing angles, plus an extra CP phase: 


= -2 
fie 


f Am\,L\ ( AmoiL 
sm( 20 i 3 ) sm sm ( j sm i 


\ 4:E, 


31“ 1 / 01“ \ T T 

+ J^NP — 1123 


'Am? L 

- cos 6*13 cos 6*23 sin( 26 *i 2 ) sin ( 


4P,/ 


sin(/jvp), 


(23) 


with I 123 = -5cp = 012 - 013 + 023 and Inp = 0i2 - Arg(Q; 2 i). 

Notice that the conversion probability depends on just two phases, the standard one, 
-6123 = and another phase describing the new physics, Inp- This new phase contains the 
information of the imaginary part of 0 : 21 , that is, the overall effect of all the additional phases 


associated with the heavy states. Notice that, besides the standard CP term in Eq. (23), two 


new CP phase-dependent terms appear; the first involves the difference between standard 
and non standard phase: /123 — Inpi while the second one depends only on Inp- One sees in 


Eq. (23) that the hrst term is proportional to sin 0 i 3 , while the second one depends on the 


solar mass difference Am^i and, therefore, both terms should be small. In order to illustrate 
their impact upon current neutrino data analysis, we show in Fig. how this new phase 
parameter influences the conversion probability. In this hgure we compare the standard 
three neutrino probability (with a “best-£t” phase 5 = —/123 = 37 r/ 2 ), with the case of an 
additional neutral heavy lepton with overall contribution given by an = 1, 0:22 = 0.9997, 
|q; 2 i| = 0.0264, and for the particular new physics phase parameter of either 7 r /2 or 37 r /2 
(left panel) or 0 , 7 r (right panel). One sees that the effect of the additional phase in future 
oscillation appearance experiments could be sizeable and, depending on the specihc value of 
this new phase, the survival probability could either increase or decrease. 

For the sake of completeness, we also give the expression for the survival probability 




cos 


j>i 


I 2E 


(24) 
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FIG. 3: Conversion probability for a fixed neutrino energy Ey = 1 GeV. The solid (black) curve 
shows the standard conversion probability, with 5 = —/123 = 37r/2. The non-unitary case is 
illustrated for an = 1, q ;22 = 0.9997, and |q! 2 i| = 0.0264. In the left panel, two values for the 
new CP phase parameter I]\fp are considered: tt/2 (dashed/magenta line) and 37r/2 (dotted/green 
line), while in the right panel we take Inp = 0 (dashed/magenta line) and tt (dotted/green line). 


P = 


(10211'^ + - 4^ |iV„f 

j>i ^ / 


Pfifj, = (|Q:2iP + (^22^ ~ \0i2lUej + H + tt22^/xiPsin^ f Lj 

j>i ^ ^ 

SO that, neglecting cubic products of 0 : 21 , sin^is, and Am^i, we will obtain 

PiM/i = Oi22P^^^ + + 2|a2i|^aL-P/!M 

with the standard oscillation formula, given by: 


(25) 

(26) 


(27) 


p3x3 ^ 


1-4 [cos^ 023 sin^ 023 - cos( 26'23) sin^ 023 sin^ sin^ f 


f AmhL 

AE 


( Arn? L\ [Am? L 

—j sin f 


-4 


( Am? L 

—) + cos^ 012 cos'^ 023 sin^ 012 


2E 


sin 


2E 

2 ( Aml^L 


2E 


K 


AE 


( 28 ) 
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FIG. 4: Correction to the standard muon neutrino survival probability for different values of the 
new CP phase parameter Inp, with the remaining parameters hxed as in Fig. 


while the extra terms in the oscillation probability are given by: 


ph ~ 


[sin 9 i 3 sin 623 cos( 26 ' 23 ) cos(Ii 23 - Inp)] sin^ 


/ AmoiL\ 


+ 2 



(29) 


Pll « 1 - 2 sin" sin" • (30) 

As for the conversion probability, —)■ Ve), we also compute the muon neutrino 

survival probability and show its behaviour in Fig. As one can see, this disappearance 
channel is also sensitive to the new CP phase. The computations were performed for the same 
parameter values used in the previous figure, that is, an = 1, 022 = 0.9997, |a 2 i| = 0.0264, 
and an overall phase of either 7 r /2 or 37 r /2 as well as 0 or tt. The Standard Model phase was 
hxed to be 5 = —/123 = 37 r/ 2 ). 


We now turn our attention to oscillations of electron neutrinos or anti-neutrinos relevant, 
say, for the description of solar neutrino experiments, as well as terrestrial experiments using 
reactors or radioactive sources. The electron (anti) neutrino survival probability (in vacuum) 
is given by the following expression: 

Fee = \Ne^?\Nei\^ + COS ( 1 , (31) 

i j>i ^ ' 
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( 32 ) 


and, using Eq. (|^, it is easy to see that N^i = anUei which leads to the expression 


Pee ^11 


|f/e*n^e*P + X;2|f^einf/e*PcOS ( 


/ Am?, 




j>i 


2E 


This transforms, in a straightforward way, to the equation 


Pee = ati [l - cos"^ ^13 sin2(26'i2) sin 2 (Ai 2 ) - sin2(26'i3) sin^(Ai 3 )] , (33) 

A m ?. 

with Aij = -^^L. Notice that in this case, the effect of a neutral heavy lepton will be an 
overall factor that accounts for the violation of unitarity: unlikely to produce visible 

effects in oscillations of, say, reactor neutrinos, given the strong universality restrictions 
derived in Fig. 

For completeness we mention that, should the extra neutrino states be light enough to 
take part in oscillations, they could potentially play a role [SSI [SS] in the anomalies reported 
by the MiniBooNE collaboration [22] or the reactor neutrino experiments 157|. We will not 
consider this possibility here. 


VI. BOUNDS FROM NEUTRINO OSCILLATION EXPERIMENTS 

From the previous formulas for the oscillation probabilities one sees that, even at zero 
distance, the survival and conversion probabilities differ from one and zero, respectively. 
This is a well-known behaviour and it is a consequence of the effective non-unitarity of the 
3x3 leptonic mixing matrix [SB]- We can express these probabilities, for the zero distance 
case, as 

Fee = < = [(iViVt)n]' = [1 - {SS^)n? , 

P,, = (|« 2 ir + « 22 )^ = = [1 - , (34) 

P,e = «?il«2ir = [{NN%,]^ = [(FFt) 2 l]l 

In order to make a quick estimate of the constraints on the new parameters, we write these 
expressions in a different way, in order to compare them with the corresponding expressions 
for a light sterile neutrino in the limit of Am‘fjL/{AE) 3> 1 ((sin^(Am?^L/(4F))) = 1/2). 
The result for our case can be expressed in an analogous way as in the case of extra light 
neutrinos [59] : 

^ee = 1- ^ [sin2(20ee)]gff, 

p,, = l-^[sin2(2M]efr> (35) 

F^e = ^ [sin 2 ( 20 ^,)]^jj. 
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with 


[sin2(20ee)],g = 2(1 

[sin^ ( 26 ^^)]^^ = 2[1 - (|a2iP + ^ 22 )^], (36) 

[sin^ (20^e)]gg = 2al^\a2i\'^. 

We can compare these expressions with the current constraints on light sterile neutrinos in 
order to get the following 3 ct limits [SB] 

[sin 2 ( 20 ee)],ff < 0 . 2 , 

[sin2(20^^)]^g < 0.06, (37) 

[sin2(20^e)],g < 1 X 10-^ 

However, we prefer to use the bound from the NOMAD experiment [BU], since it is the 
most reliable constraint on the zero-distance effect (neutrino non-orthonormality due to 
heavy neutrino admixture) from neutrino oscillations. Translated into the parametrization 
under discussion, this constraint takes the form 


ah\a2i\^ < 0.0007 (90%C.L.) 


(38) 


If we combine this limit with those coming from universality at Eqs. (12) and (15), the 
following 90% C.L. bounds (1 d.o.f.) are obtained 


0.989, ai2> 0.999, |a2ip < 0.0007. 


(39) 


VII. COMPILING CURRENT NHL CONSTRAINTS 

Non-standard features such as unitarity violation in neutrino mixing could signal new 
physics responsible for neutrino mass. For example, they could shed light upon the prop¬ 
erties of neutral heavy leptons such as right-handed neutrinos, which are the messengers 
of neutrino mass generation postulated in seesaw schemes. In many such schemes the 
smallness of neutrino masses severely restricts the magnitudes of the expected NHL sig¬ 
natures. However these limitations can be circumvented within a broad class of low-scale 
seesaw realizations [T^l - EU] . For this reason in this section we will present a compilation of 
model-independent NHL limits, which do not require them to play the role of neutrino mass 
messenger in any particular seesaw scheme. Results of this section are not original, but they 
are included for logical completeness. 

Isosinglet neutrinos have been searched for in a variety of experiments. For example, if 
they are very light they may be emitted in weak decays of pious and kaons. Heavier ones. 
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but lighter than the Z boson, would have been copiously produced in the hrst phase of the 
LEP experiment should the coupling be appreciable [25l |26]. Searches have been negative, 
including those performed at the higher, second phase energies |27j . 



m (GeV) 


FIG. 5: 


Bounds on the component of a heavy isosinglet lepton of mass in the electron neutrino. 


A summary of constraints for the direct production of neutral heavy leptons is shown in 
Figs. 00 and In most cases, experiments have looked for a resonance in a given energy 
window, for a given mixing of the additional state, described in this case by the submatrix 
S of Eq. Q. Although the constraints for the mixing in these cases are stronger, in most of 
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FIG. 6: Bounds on the component of a heavy isosinglet lepton of mass rrij in the muon neutrino. 
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FIG. 7: Bounds on the component of a heavy isosinglet lepton of mass rrij in the tau neutrino. 

the cases they rely upon extra assumptions on how the heavy neutrino should decay. 

In particular, in Fig. we summarize the constraints on \Sej\^ for a mass range from 
10“^ to 10^ GeV coming from the experiments TRIUMF [HTl (denoted as vr — )• ei/ and 
R —)■ ez/ in the plot), PS191 [63], NAS jH], CHARM [65], Belle [66], the LEP experiments 
DELPHI [27], L3 [HZ], LEP2 [6B], and the recent LHC results from ATLAS [HSIZD]. Future 
experimental proposals, such as DUNE mi and ILC, expect to improve these constraints mi 
In Fig. [^ we show the corresponding constraints for the case of the mixing of a neutral 
heavy lepton with a muon neutrino. In this case we show the experimental results coming 
again from PS 191, NAS, and Belle, from the LEP experiments L3, DELPHI, and from the 
LHC experiment ATLAS; we also show the bounds coming from KEK [73l [71] (denoted as 
K ^ fiu in the plot), CHARM H [75], FMMF [76], BEBC [77], NuTeV [78], E949 [79], and 
from the LHC experiments CMS [HO] and LHCb [HI]. Finally, for the less studied case of the 
mixing of a neutral heavy lepton with a tau neutrino, the known constraints, coming from 
NOMAD [H 2 ], CHARM [HH], and DELPHI [27] are shown in Fig. ([^. 

Heavier neutrinos in the TeV range, natural in the context of low-scale seesaw, can also be 
searched for at the LHC. However, within the standard SU(3)c ® SU(2 )l ® U( 1 )y model such 
heavy, mainly isosinglet, neutrinos would be produced only through small mixing effects. 
Indeed, it can be seen from Figs. and that restrictions are rather weak. In contrast, 
this limitation can be avoided in extended electroweak models. In such case a production 
portal involving extra kinematically accessible gauge bosons, such as those associated with 
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left-right symmetric models, can give rise to signatures at high energies, such as processes 
with lepton flavour violation (HUES]. 

A. Neutrinoless double beta decay 

If neutrinos have Majorana nature, as expected on theoretical grounds, neutrinoless dou¬ 
ble beta decay is expected to occur at some level [30] • We start our discussion by reminding 
the dehnition of the effective Majorana neutrino mass |86j . 

(m) = I (40) 

3 

where the sum runs only for the light neutrinos coupling to the electron and the W-boson. 

From Eq. ([^ one sees that, in the presence of the heavy neutrinos, the three light SM 
neutrino charged current couplings will be modihed to = cxnt/ei, with i=l,2,3, and 

their contribution to neutrinoless double beta decay will change correspondingly. 



FIG. 8: Sensitivity of neutrinoless double beta decay to isosinglet mass rrij in the electron neutrino. 

Moreover, the heavy states will induce also a short-range or contact contribution to neu¬ 
trinoless double beta decay involving the exchange of the heavy Majorana neutrinos. Since 
these are SU(3)c ® SU(2 )l ® U(1)y singlets they couple only through the mixing coefficients 
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Sej- The general form of the amplitude is proportional to 


A oc 


9 2 ’ 

— rrij 


( 41 ) 


where q is the virtual neutrino momentum transfer. Clearly there are two main regimes for 
this amplitude; for we have 



(42) 

1 

*^heavy ^ 

m,- 

(43) 


while for 


This behaviour can be seen in the corresponding estimated sensitivity curve shown in Fig. 
(|^. This line is obtained for ^®Ge assuming a single massive isosinglet neutrino [HZ|- The 
change in slope takes place for masses close to the typical nuclear momentum, around 
100 — 200 MeV. Both light and heavy contributions must be folded in with the appropriate 
nuclear matrix elements [SB] whose uncertainties are still large. As a result it is not possible 
to probe the indirect NHL effect upon the light neutrino contribution to the effective mass 


in Eq. (40) which amounts to a multiplicative factor in the amplitude, a difference well 


below current sensitivities. Notice that, in contrast to bounds discussed in Figs. |^[^and[^ 
the restriction from the neutrinoless double beta decay in Fig. holds only if neutrinos have 
Majorana nature. 


B. Charged lepton flavour violation 

Virtual exchange of NHLs would also induce charged lepton flavour violation processes 
both at low energies [29] as well as in the high energies provided by accelerator experi¬ 
ments [33] . However rates would depend on additional flavor parameters and upon details 
on the seesaw mechanism providing masses to neutrinos. The possibility of probing it at 
hadronic colliders such as the LHC may be realistic in low-scale seesaw models with ad¬ 
ditional TeV scale gauge bosons beyond those of the SM gauge structure and with lighter 
NHLs [HU lESl ISni ED] . However we do not consider this possibility any further here because 
the corresponding rates depend on very model-dependent assumptions. 


VIII. SUMMARY 

Simplest seesaw extensions of the Standard Model predict unitarity deviations in the 
leptonic mixing matrix describing the charged current leptonic weak interaction. This is due 
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to the admixture of heavy isosinglet neutrinos, such as “right-handed neutrinos”, which are 
the “messengers” whose exchange generates small neutrino masses. Low-scale realizations of 
such schemes suggest that such NHL may be light enough as to be accessible at high energy 
colliders such as the LHC or, indirectly, induce sizeable unitarity deviations in the “effective” 
lepton mixing matrix. In this paper we used the general symmetric parametrization of 
lepton mixing of Ref. [8] in order to derive a simple description of unitarity deviations in 
the light neutrino sector. Most experiments employ neutrinos or anti-neutrinos of the hrst 
two generations. Their description becomes especially simple in our method, Eq. ([^, as 
it involves only a subset of parameters consisting of three real effective parameters plus a 
single CP phase. We have illustrated the impact of non-unitary lepton mixing on weak decay 
processes as well as neutrino oscillations. For logical completeness we have also re-compiled 
the current model-independent constraints on heavy neutrino coupling parameters arising 
from various experiments in this notation. In short, our method will be useful in a joint 
description of NHL searches as well as upcoming precision neutrino oscillation studies, and 
will hopefully contribute to shed light on the possible seesaw origin of neutrino mass. 
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IX. APPENDIX: NEUTRINO MIXING AND HEAVY ISOSINGLETS 

As already explained, heavy gauge singlet neutrinos arise naturally in several extensions 
of the Standard Model. The general form of the mixing matrix describing their charged 
current weak interaction has been given in [8]. Here we will further develop the formalism 
so as to describe not only the couplings of the additional heavy neutrinos but also their 
effects in the light neutrino sector in a convenient but complete way, with no assumptions 
about CP conservation. Using Okubo’s notation [HT], we can construct the rotation matrix 

jjnxn 


u 


nxn 


— UJn-ln^n-2n ••• ^In ^n-2n-l ^n-Sn-1 ••• ^ln-1 ••• <^23 <^13 ^12? 


(44) 
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where each ojij {i < j) stands for the usual complex rotation matrix in the ij plane 

^ ci3 0 ^ 

a;i3 = 0 10 , (45) 

\ 0 ci3 / 

with Sij = sin^jj and Cij = cos9ij. This matrix can be expressed in general as: 

^ 0 / 3 pj^ij T T TJij^aj^jSi j (^6) 

where i < j and sfj = sin^ 9ij, rjij = sin6*jj and sin6*jj, generalizing the 


matrix in Eq. (45) as: 


/l 0 ■■■0 ■■■ 

0 1 


UJij — 


\ 




0 ■■■ Tlij 


0 


9ij ' ' ' 0 • • • Cij 


0 


1 0 
0 1 


In general, one can decompose Eq. (44) in the following way 


jjnxn _ jjn—N JJ^ 


with 


— OJm-IN ^N-2N ■ ■ ■ ^IN, 


Y jn _iV 

U ^n—ln ^n—2n • • • ^In ^n—ln—1 ^n—2n—l • • • ^ITV+l? 

SO that the matrix decomposition will be given by 

/ 


jjn—NjjN _ 


011 0 ■ ■ ■ 0 

0(21 Oi22 ■ • • 

! 0 

CtNl ■ ■ ■ Q-NN 


V 


\ ( ttN ttN ... ttN 


s 


'^11 ‘^12 '^IN 

tjN tjN : 

'^21 *^22 


tjN 

'-'ni 


ttN 

'-‘nn 


(47) 


(48) 

(49) 

(50) 


\ 


( 51 ) 
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which turns out to be very convenient. The 3x3 neutrino mixing matrix, determined 

in oscillation experiments could be unitary, or it could be just a non-unitary submatrix of 


the larger mixing matrix described in Eq. (44). Therefore, when dealing with more 


than three neutrinos, we can write as the product of two matrices: 


jj-nxn _ jj 


NP jjSM 


(52) 


where "NP" means ’’new physics” and "SM" stands for the “Standard Model” matrix, 

= (jJn-ln^n-2n ■ ■ ■ ^‘in^2n ^ln^n-2n-l ■ ■ ■ ^2n-l ^In-l ■ ■ ■ 1^34 ^2A , (53) 

JjSM = 3 (Ui 3 (Ui 2 . (54) 

The complete n x n matrix, may be written as [35] 

N S 


jjnxn _ 


1/ T 


(55) 


where N is the 3x3 matrix with the standard neutrino terms. From Eq. (52) one sees that 
N can always be parametrized as 


N = = 


^ dll 0 0 

0^21 0^22 0 

V ttsi tt32 tt33 / 




3x3 


(56) 


where the zero triangle submatrix characterizes this decomposition. It is useful to see how 
the components of this matrix can be found. First notice that oOijOOki commutes when 


i ^ k,l and j ^ /c, /; therefore, Eq. (53) can be rewritten as 

jjNP _ 

U ^n—2n • • • ^An ^n—2n—l ■ ■ ■ ^An—1 ■ ■ ■ ^45 ^ 

^2n ^In ^3n-l <^2 n-1 n-1 ■■■ ‘^SA 1^2A^1A- 


(57) 


Clearly, the first line of this equation has no influence in the submatrices N and S. On the 
other hand, the second line of the above equation is a set of products of the form LJ 3 jLJ 2 j‘^ij, 
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each of them having the form: 


— UJ2,jOJ2jOJij 


Cij 

0 

0 


0 


Vij 

0 


V2jVij 

C2j 

0 


0 


V2jClj 

0 


V3jC2jVlj 

W2j 

C3j 


0 

V3jC2jClj 

0 


0 

0 

0 


I 


0 

0 


C3jC2jVlj 

C3jf]2j 

V3j 


0 

C3jC2jCij 

0 


0 

0 

0 


0 


0 

I 

/ 

a{i 0 

0 ; 

0 

a{j 

0 

\ 




*^21 *^22 

0 : 

0 

c^ij 

0 





<^31 *^32 

^33 ■ 

0 

a3j 

0 





0 

0 

0 ; 

I 

0 

0 





«2i 

“i3 • 

0 

• 

jj 

0 





0 

0 

0 ; 

0 

0 

I 

) 





(58) 


We can see that the expression for depends only on products of the type 

Q,nQ,n-i... After performing the multiplication one notes that the diagonal entries 

of the matrix are in general given by 



n 

n —1 

n —2 

4 




ttll 

= «ii 

«ii 

«ii ■■ 

■ «ii = 

Cl n 

CIn—iCln— 2 •* 

• Cl 4 


n 

n —1 

n —2 

4 




tt 22 

— ^22 

•^22 

0^22 

■ ^22 — 

C2 n 

C 2 n-lC 2 n -2 • • 

• C24 


n 

n —1 

n —2 

4 




0^33 

— <^33 

^33 

“33 ■ 

■ ■ ^33 — 

C3 n 

CSn—1C3n—2 • 

• • C34 


while the off-diagonal entries ck^ are given as: 


«21 

= 

•^21 

«ir' 

4 

■■■ail 

+ 

n n—1 

^22 ^21 

4 

■■■ail 

+ ■ 

■■■ + 

^22 

n—1 n—2 

^22 ^22 

4 

C^21 5 

«32 

= 

•^32 

«22“^ 

4 

CI 22 

+ 

^ n ^ n—1 
•^33 *^32 

4 

CI 22 

+ ■ 

■■■ + 

•^33 

^ n—1 ^ n—2 

“33 “33 ■ ■ ■ 

4 

*^32 7 

0^31 

= 

•^31 

till 

4 

■■■ «ii 

+ 

^ n ^ n—1 
•^33 '^31 

4 

■■■ail 

+ ■ 

■■■ + 

“33 

^ n—1 ^ n—2 

“33 “33 ■ ■ ■ 

4 

^^31 


+ 

•^32 

( ^ n—1 
( *^21 

ail 


ail *^22 

l^n-2 

•^21 


4 1 

ail H- 

' fo CI 22 *^22 

4 

*^21 


+ 

+ 

+ 


n 

33 


(«2r^ 

ail 

4 

■■ail 

+ ■■■ + 

n—2 n—3 
^22 ^22 

■ ■ ■ a: 

n 

33 


«32-^( 

a^f ^ a 

n —4 

11 

4 1 

■ ■ “11 + ■ 

■ ■ ■ + a22 ^ 

^ n —4 
*^22 

n 

^ n—1 

^ n—2 


^4 




33 

*^33 

“33 

^32 

0^21 ) 





) + 


a 


21 
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or, more explicitly, 


<^21 = C 2 n C 2 n-l • • • C 25 ^24^14 + C 2 n ■ ■ ■ C 26 ^25^15 C 14 + ... + ?72n%n Cln -1 Ci „_2 ... C 14 , 

<^32 = CSn Csn-l . . . C 35 ?734?724 + CSn ■ ■ ■ C 36 ^35^25 C 24 + • • • + %n^2n C2n-1 C2n-2 • • • C 24 , 

<^31 = CSn C3n-1 . . . C 35 ?734C24^14 + Csn . . . C 36 ^35C25^15 C 14 + . . . + r]^nC 2 n Vln Cln -1 Cln -2 ■ ■ ■ C 14 

+ C 3 „ C 3 „_i . . . C 35 ? 735 r 725 ^ 24%4 + C 3 „ . . . C 3 g ^736^26C25 ^24^4 

+ • • • + V3nV2nV2n-lVln-lCln-2 . . . C 14 . (60) 


With these formulas, and the known expression for we already have the explicit 

description of Eq. (56) for any number of extra neutrino states. Before concluding this 
appendix, we would like to remark that the position of the three off-diagonal zeros in 
was chosen to conveniently make the matrix lower triangular. This simplihes the form of the 
non-unitary lepton mixing matrix describing most situations of phenomenological interest, 
involving solar, atmospheric, reactor and accelerator neutrinos. By choosing alternative 
factor-orderings, one can have different parameterizations, with the zeros located at different 
off-diagonal entries. 


Application to 3 -|- 1 seesaw scheme 


We will conclude this appendix by showing the expressions for aij in the case of one and 
three additional neutrinos. For the case of just one additional neutrino, the mixing matrix 
is given by 

^4x4 ^ I ^3x3 *S'3xl \ 

\ Tix 3 Vlxl j 

The corresponding expressions for the parameters aij will be given by 



«11 = Ci4 , 

0^22 = C24 , 

«33 = C 34 , 

«2i = h24 Via ; 

«32 = VSA V2A , 

0^31 = VsA C24 Via ■ 


(62) 
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Application to 3 + 3 seesaw scheme 


In this case the full mixing matrix will have the following structure 


jj6x6 

with the a parameters given by 


-^3x3 

^3x3 


>S'3x3 

h"3x3 


«11 = C16C15C14, 

«22 = C26 C25 C24 , 

«33 = C36 C35 C34 , 

«21 = V26 VI6 Cl5 Cl4 + C26 Vl5 Cl4 + C26 C25 ?724 Vl4 , 

«32 = C36 C35 7^34 f/24 + C36 C35 f]25 C24 + ^36 V26 C25 C24 , 

0^31 = C36 C35 C34 7734 C24 t/l4 + C36 7^35 C24 t/l5 Cl4 + V36 C26 hl6 C15 C14 

+ C36 7^35 7^25 ^24 ^14 + ^36 V26 C25 ^24 i?14 + ^36 V26 V25 hl5 C14 • 


( 63 ) 
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